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Quasiparticle states bound to the boundary of anisotropically paired superconductors, their 
contributions to the density of states and to the conductance of NIS junctions are studied both 
analytically and numerically. For smooth surfaces and real order parameter we find some general 
results for the bound state energies. In particular, we show that under fairly general conditions 
quasiparticle states with nonzero energies exist for momentum directions within a narrow region 
around the surface normal. The energy dispersion of the bound states always has an extremum 
for the direction along the normal. Along with the zero-bias anomaly due to midgap states, we 
find, for quasi two-dimensional materials, additional low-temperature peaks in the conductance of 
NIS junctions for voltages determined by the extrema of the bound state energies. The influence 
of interface roughness on the conductance is investigated within the framework of Ovchinnikov's 
model. We show that nonzero-bias peaks at low temperatures may give information on the order 
parameter in the bulk, even though it is suppressed at the surface. 

PACS numbers: 74.50. +r, 74.80. F, 74.72.-h 



I. INTRODUCTION 

It is now widely recognized, that the anisotropy ratio of the order parameter on the Fermi surface for many of 
the high- 2^, superconductors is of the order of unity. By contrast, for conventional low-temperature superconductors 
the anisotropy ratio is known to be much smaller, so that they are almost isotropic s-wave superconductors. This 
fact itself, essentially irrespective of the specific type of the superconducting pairing, results in many new important 
consequences for the properties of HTSC, compared to isotropic s-wave superconductors. In particular, a highly 
anisotropic order parameter turns out to be quite sensitive to any kind of inhomogeneities in the material, including 
nonmagnetic impurities as well as interfaces. 

For most surfaces the bulk behavior of an anisotropic order parameter can differ essentially from its behavior at 
the wall. Even smooth specularly reflecting walls are, in general, pair breaking, which results in a suppression of 
the order parameter near the boundary. This concerns both anisotropic superconductors and superfluids.EJiD At the 
same time, quasiparticle states bound to the boundary appear^du£-.to this suppression or/and due to the sign change 
of the anisotropic order parameter with momentum direction.oaiij Furthermore, in the region where the bulk order 
parameter is essentially suppressed, subdominant pairing channels with different symmetries may come into play and 
become stable near the wall.0u 

Several important experimental methods used for studying the anisotropic structure of the order parameter, in 
particular tunneling measurements, are fairly sensitive to the superconducting properties close to the surface of the 
sample. The local quasiparticle spectrum at the surface which differs, in general, from the bulk density of states is of 
crucial importance for the I-V characteristics of tunnel junctions. Apart from the fact, that the gap anisotropy washes 
out those peaks in the tunneling current which would occur in the isotropic case, the presence of quasiparticle bound 
states or/and the order parameter of subdominant symmetry channels _LLear the surface results in new characteristic 
features of the I-V curves in the vicinity of some specific points. One of the most striking features of NIS 
junctions with^-wave superconductors is the zero-bias anomaly of the conductance which is associated with surface 
midgap states.EJ These states localized at the surface are a consequence of the sign change of the order parameter 
on the FefPH surface. A further consequence of those states is a low-temperature anomaly in the Josephson critical 
current .E3ilJ By contrast, some other important features of the Josephson effect do not depend on the quasiparticle 
spectrum at the surface, having a more universal character. Actually, they cLeeeBd only on the phase difference of the 
order parameters in the bulk superconductors on both sides of the junctionJialill 

Below we develop an analytical investigation, including numerical calculations, of quasiparticle states bound to the 
boundary of anisotropically paired superconductors, their contributions to the density of states and to the conductance 



1 



of NIS junctions with-weak transparency. Our approach to the problems is based on the quasiclassical formulation 
of superconductivity.EaE^I Both smooth specularly reflecting and diffusive interfaces are considered. The spatially 
dependent order parameter is supposed to be real. The analytical study of quasiparticle states with nonzero energies 
bound to a smooth wall is carried out for the first time, and some general results are obtained. In particular, we 
show that the energy of a bound state is closely connected with the phase of the complex residue of the anomalous 
propagator. Quasiparticle bound states with nonzero energies turn out to exist under fairly general conditions for 
momentum directions within a narrow region around the surface normal for which the energy dispersion always has 
an extremum. 

While the presence of a zero-bias anomaly in the conductance does not depend upon details of the Fermi surface 
but only on the appearance of midgap states, low-temperature anomalies in the conductance due to nonzero-energy 
levels near the surface are absent for three dimensional systems. For quasi two-dimensional materials we find, that 
quasiparticle bound states with nonzero energies lead, at least at low enough temperatures, to peaks in the conductance 
for voltages determined by the extrema of the bound state energies. It is shown below, that for a smooth barrier plane 
the magnitudes of those peaks are proportional to the inverse square root of the temperature (G(T — > 0) cx 1/VT) in 
contrast to zero-bias anomalies for which G(T — > 0) cx 1/T. 

In realistic systems, the peaks in the quasiparticle spectrum will be broadened. The shape and the height of 
these peaks are governed by the physical parameters of the junction like interface roughness (e.g., lattice appjtant 
mismatch) and transparency. We investigate the influence of interfacfi-caughness on the local density of states&C-a and 
on the conductance within the framework of Ovchinnikov's model.BoO In all our calculations the order parameter 
is determined self-consistently. In contrast to midgap states, nonself-consistent calculations, disregarding surface pair 
breaking, do not permit to find and describe quasiparticle bound states with nonzero energy and the corresponding 
peaks in the conductance at nonzero-bias voltage. 



II. QUASICLASSICAL THEORY 
A. Basic Equations 

In general, strongly anisotropic superconductors react sensitively to surfaces and interfaces. For instance, the 
superconducting order parameter near the surface is reduced due to pair-breaking effects. The quasiclassical theory of 
superconductivity is the most efficient theory for investigating effects of surfaces and interfaces. We use Eilenberger's 
equations for the quasiclassical retarded propagator g R which for a clean singlet anisotropically paired superconductor 
reduce to the following 2x2 matrix form:E3 

ef 3 - A(p f ,R),g R (p f ,R;ej\ + iv/-V K g R {p f , R, e) = , (I) 
[g R (p f ,R;s)] 2 = -7r 2 l . (2) 

Here, e is the quasiparticle energy, pf the momentum on the Fermi surface, Vf the Fermi velocity, and A the order 
parameter matrix. Following standard notations we use a "hat" to indicate matrices in Nambu space. A convenient 
basis set of Nambu matrices is the unit matrix I and the three Pauli matrices fx, f%, T3. The propagator g (here and 
in the following we drop the superscript R for simplicity) and the order parameter matrix A have the form 

9 ^ and A=( ° . (3) 



/+ -g J " \ -A* 

The gap function A(p/, R) is related to the anomalous Green function / and has to be determined self-consistcntly 
(see Ref. Q) . The diagonal part g of the full matrix propagator g carries information on the superconducting density 
of states (DOS), 

v(p f ,R;e) = Jm[g(p f ,R]s)] . (4) 

This quantity contains all the information we need for calculating the tunneling current across an NIS interface (see 
Eqs. pi)). ' 
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B. Interface Model 



The boundary conditions ai-an ideal interface are given by Zaitsev's relations (see Refs. pq,0), which reduce in the 



limit of zero-transparency tol 



g(pf) = g(p f ) , (5) 



where the propagators are taken at the metal-insulator boundary. Equation (|5|) connects the propagator of an incoming 
quasiparticle with Fermi momentum p j and the propagator of the reflected quasiparticle with Fermi momentum p ^ at 

the interface. For specular reflection the momentum parallel to the interface is conserved, i.e., pi = pi. For a complete 
determination of the quasiclassical propagator we have to take into account that deep inside the superconductor the 
propagator approaches its bulk value. 

In realistic systems the reflection of quasiparticles fro^sttrlaces and interfaces is expected to be diffuse, at least 
to some degree, due to imperfections at the boundary.BEj'Elo For isotropic superconductors, elastic scattering of 
quasiparticles on walls has no effect on the density of states near the interface, and thus, on the tunneling current. 
The reason for this is that the order parameter for the reflected quasiparticles is the same for all directions. In 
anisotropic superconductors the value of the order parameter is no longer independent on direction. A quasiparticle 
on an outgoing trajectory will see an order parameter which is, in general, different from the one on its incoming 
trajectory. This leads to pair breaking even for a specularly reflecting wall, and, as a consequence, to a suppression 
of the order parameter near the boundary. For diffuse scattering of quasiparticles at the interface all outgoing 
trajectories are mixed. As a result, sharp structures in the density of states (e.g., bound states) which would occur 
for ideal interfaces are smeared out near rough interfaces. Thus, the effect of surfaces and interfaces is similar to 
the effect of static impurities which leave the order parameter of isotropic superconductors unaffected ("Anderson's 
theorem"), but have destructive influence in anisotropic systems where superconductivity can even disappear when 
the impurity concentration exceeds a critical value. 

We model a rough interface by coating an ideal interface (transparency D) on both sides by thin layers of strongly 
disordered metals which differ from the respective bulk materials only by their very short mean free paths (see Fig. 
1). The roughness of the interface is measured by the phenomenological parameter p which is essentially the ratio of 
the thickness of the layers 5 and the mean free path £, i.e., p — 5/1. The case p = corresponds to an ideal interface, 
whereas p — oo describes a diffuse interface. The transport of quasiparticles across the disordered layers is determined 
by the following transport equationnJ 



The superscripts / and r stand for the left and right sides of the interface. The dimensionless velocity v f' r is the 



2n (9 l -%>9 l - r \ + iv^ r d x g l ' r = . (6) 
left and rig 

perpendicular to the interface component of the quasiparticle velocity in the bulk material normalized by an averaged 

Fermi velocity, v l,r = vj' 1 */ ^/(l v /' T 'l 2 )p / • ^ e nrs t term in Eq. (^j) is Ovchinnikov's anisotropic scattering self-energy 

for incoming and outgoing quasiparticles .El At the ideal interface which separates the two layers the quasiclassical 
propagators are discontinuous and the jump is given by Zaitsev's relations.Ba The normal state resistance Rm of the 
interface turns out to be independent of the roughness and is given by 

v', ± >0 

where A is the area of the interface. A detailed description of this interface model can be found in Ref. [lj. 



III. QUASIPARTICLE BOUND STATES AT SPECULARLY REFLECTING WALLS 

In this section, we study in detail the appearance of quasiparticle bound states near specularly reflecting impene- 
trable surfaces. Bound states are localized within a few (or even many; see below) coherence lengths near the surface, 
and their energies lie below the bulk value of the (momentum dependent) gap. Their existence manifests in <5-peaks 
in the momentum resolved density of states (^|) below the continuous spectrum. There are two types of bound states 
of different origin: a) zero-energy bound states or midgap stateso which are a consequence of the change of sign of 
the order parameter along a quasiparticle trajectory which is reflected from the wall, and b) nonzero-energy bound 
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states which emerge due to the depletion of the order parameter associated with the interface.!!™] While zero-energy 
bound states are a robust phenomenon of quite general origin,E3 less information is available for nonzero-energy bound 
states. 



A. Analysis of Quasiparticle Bound States at Specular Walls 

Let the quasiclassical propagator g have a pole at e = £ B (pf), where £ B (pf) is the energy of a bound state with 
Fermi momentum pf (note that E B (pf) — £ B (p ^))- It is convenient to introduce the residue of the propagator g by 

g{p f ,R;e B (pf))= Hm [(e - e B (jpf))g{p h R;e)] , (8) 

which is finite and satisfies the same transport equation as g, but with a modified normalization condition, 

[l(p f ,R;s B (p f ))] 2 =0 . (9) 

The boundary conditions for g coincide with (||). But, in contrast to g the new propagator g vanishes in the bulk of 
the superconductor. Thus, the quantity g describes just a quasiparticle state bound to the boundary plane. 

For the sake of dcfmitcncss, the superconductor is supposed to occupy the half-space x > with an impenetrable 
boundary. Furthermore, we assume that the gap function A can be chosen real. Then one obtains from Eilenberger's 
equations: 

f + {p f ,x;£ B (p f )) = -f*(pf,x;s B (p f )) , 
g(p f ,x;e B (pf))=g*(p f ,x;e B {pf)) = \f(pf,x;e B {pf))\ . (10) 

The last identity is a consequence of Eq. (^) . We introduce the phase of the complex anomalous Green function by 

f(pf,x-,s B (pf)) = -g(p f ,x;e B (pf))exp(itp(pf,x)) , (11) 

and obtain from (|l]), (||) 



g(pf,x;£ B (j3f)) = g(pf,0;e B {pf))exp I J A (Pf, sin (<p(p f ,x)) dx J , 



(12) 



together with the following equation for the phase, 



V f 

- -^fd x ip(p f ,x) +e B (p f ) - A(p f ,x) cos (ip(p f ,x)) = , (13) 

and the boundary and asymptotic conditions, 

<p(pf, 0) = (p(p , 0) and v/ lX Aoo(p/) sin (Voq(p/)) > . (14) 

If ip(pf,x) is a solution of Eqs. (|l3|), ( |l4"|) with the energy £ B (pf), then ir — ip(pf,x) is a solution with the energy 
- £ b(p/). _ 

Since d x (p(pf,x) vanishes in the bulk (x — > oo), we immediately get from Eqs. ( |l3| ) and (14h the relation between 
the bound state energy £ B (jJf) and the phase in the bulk of the superconductor, 

?B{Pf) = Aoo^COS^ooCP/)) = A 00 (p / )cos UpcvipM . (15) 

As a consequence, bound states might exist for a given momentum direction only below the band edges for the 
momenta Pf and p , i.e., for 

\s B (p f )\ ^rmnjlA^p/^IA^ppi} . (16) 

Under certain conditions it is possible to find explicit expressions for the bound state energy e b (p/) and the phase 
(Pco(Pf) in terms of the spatial dependence of the order parameter near the surface. For this purpose we transform 
Eq. (O) into the corresponding integral equation, 
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(p(p f ,x) =(p(p f ,0) + 2A °°( p /) / ( COS^oo (p f )) - 

v 

Furthermore, we suppose that the phase (p(pf, x) deviates only weakly from its value at the surface everywhere within 
the superconductor (\(p(pf,x) — (p(pf,0)\ <C 1)- Then, as a first approximation, one can substitute on the right hand 
side of Eq. @, cos (p(pf,x)) -> cos (i/?^ (p/)) - sin (^oo (p f )) (p [x,pf) - p^ (p/)), which leads to 

A (pi) 

Voo (p/) = V {Pf, 0) + — cos (p/)) A(p f , Voo (p f )) , (18) 



*niax 



where 



oc 



V 



A(p f ,<p co (p f )) = ^m /expf--^-sin(^ 00 (p / )) j A{p f ,x)dx\ (l - A ^ P f^ )dx. (19) 



\ 



Aoo(Pf) 



Here, we introduced the quantity A max = max {| Aoo(p/)|}. 

Sf 

The condition of weak deviations of the phase from its value on the surface is transformed with the help of Eqs. 
(H), © into the form 

|A T iPf)l 1 cos ( yoo (p f )) ] \A(p f , Voo (p f ))\ « 1 . (20) 

"max 

Thus, the approximation can be justified, for example, for momentum directions for which the order parameter is only 
slightly suppressed at the surface (when \A(pf, ip^ (p/))| is small), or in the vicinity of nodes of the order parameter 
(when lAppfpf )\/A max is small ), or for a small enough | cos (poc (p/)) |. 

If condition (go|) is valid both for incoming and for outgoing momenta, we get from Eqs. (|l4|), (|l5|), ( fL8| ) the following 
equation for the bound state energy 



£fl(P/) 



2 , X 



COS 



A(p f ,i Poo (p f ))-A(p f , iPoo (2 f )) ) . (21) 



It is easy to see that the condition A 00 (pf)A 00 (p ) < implies the existence of a midgap state at esiPf) — 0. We 
should note that this solution does not depend on the validity of our expansion ( |l8| ) for the phase p^ (p/)- If the 
order parameter changes sign along a trajectory, then p(pf) — (tt/2) sgn (v/ j2 ,A 00 (p/)) is a solution of Eqs. (|l3|), (14) 
for es(p/) = 0. The left-hand side of Eq. (|2(]) is equal to zero for this solution. Another consequence of Eq. (21) 
is the absence of nonzero-energy states near the edge of the continuous spectrum (that is close to the right-hand 
side of ([l6|)) provided A 00 (p/)A 00 (p^.) < and \A(pf,cp 00 (pf))A OQ (pf)\/ A max <C 1. However, bound states may 

appear near the band edge if the parameter \A(pf,p oc (pf))A OQ (pf)\/ A max is not small (see Fig. 4c). In this case, 
our approximation ( |l8| ) for the phase p oc (p/) breaks down. 

Let us now consider the opposite case, A 00 (p/)A 00 (p^.) > 0. If the condition \A aD (pf)A(pf,p 00 (p/))|/ A max <C 1 
is satisfied for a pair of incoming and outgoing momenta, we see from Eqs. (|l^) and ( pl| ) that bound states occur only 
for trajectories with jA^py) — A 00 (p^.)| <C |Aoo(p/)|. The energies of these bound states lie close to the continuum. 
These bound states arise for momentum directions close to the normal or almost parallel to the boundary plane 
where A oa (pf) ps Aoo(p^). The latter region of momentum directions, however, is not so interesting because of its 
comparatively small contribution to the tunneling current. For bound state energies near the edge of the continuum 
(es(p/) « Aoo(p/)) we find from Eq. ( p"5| ) | sin (^^(p/)) | <C 1. Therefore, we can approximate Eq. ( |l^ ) by 



^^'•^^TA 1 - (22) 
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and obtain the following expression for bound states energies near the band edge: 

4(Pf) A 00 ( P/ )A 00 (p ) / A 0o (p / )A oo (p ) , , 2 \ 

A^pU / A 00 (p / )A 00 (p / ) ! * 




A oo(p/) / 4A 



max 



MPf)-Mp f )) J 



(23) 



It is of particular interest to consider the momentum direction along the normal to the boundary plane p^ ||n, where 
the additional relation p™ = — p^ holds. Using A(p/) = A(— p/) we get ^{plj,x) = —Lp(p r ^.,x). Together with the 
boundary condition ( |14|) fo r y> at the wall we find </?(pj,0) = 0. Then one has |Voo(P/) — ( p(p , }t^)\ = l<Poo(P/)l <C 1, 
and obtains from Eqs. (|15]) , ([l8|) the following expression for the bound state energy: 

\e B (p])\ = lAoofc,?)! cos f^M Ao (p«)) « f 1 - ^HaKp]) ) lA^p])] . (24) 

\ L^-max J y LL ^max ) 

The asymptotic condition for 95 is satisfied for this state for positive values of Ao(pJ), hence 

A(p?,:-, , 

1 - a In )dx > . (25) 




This inequality is fulfilled, in general, for the order parameter with the amplitude near the wall being less than its 
bulk value. This is true, in particular, for a one-component order parameter. Condition ( |25| ) might be violated only if 
order parameter components with different symmetries appear in the vicinity of the surface due to nonzero coupling 

constants in subdominant symmetry channels.BH In the absence of this effect, the function A i^p^x^j is always a 

monotonously increasing function of the distance from the boundary, which ensures the validity of Eq. (|25|). Thus, 
we have shown in Eq. ( |24| ) that bound states with nonzero energy do exist for Pf\\n for any (even tiny) suppression 
of the order parameter at the surface. If one interprets these states as bound states in the effective potential well 
formed by the order parameter at the boundary, one can say that for p/||ri quasiparticle bound states appear in any 
shallow potential well. 

One can find the dependence of the bound state energy on the momentum direction in the vicinity of perpen- 
dicular incidence, i.e., for pf ps p^. Assuming that the boundary is a symmetry plane of the Fermi surface of the 
superconductor, one may write 

Pf = (Pf-n)n + 5p f = p] + 5p f , p f = 6p f - p] . (26) 

Furthermore, we expand all the quantities on the right-hand side of Eq. ((2^) in powers of small deviations 5pf. In 
particular, 

A oc (p / ) = A oc (p' / l ) + (S Pr V Pf ) A oc (p^) + ~ (6 Pr X7 Pf ) 2 Aoo(p7) + . . . . (27) 
By keeping only terms up to second order in Spf, expression ( p3| ) for the bound state energy can be reduced to 

max max \ max I max 



f> 

'—^max 

2 

/ 1 A 4 frfl 



(Sp r V Pf )A^ p] )) / 1 A^(p?) 

A lo(P]) [AUp]) ^max ° {Pf) ) + 



+ (^^/a^_ 2 a^) y (28) 

The linear terms in <5p/ cancel each other, and hence the bound state energy has a local extremum for the direction 
normal to the boundary. 



G 



It follows from condition ([l6]) , that bound states occur only in a narrow region of the Fermi surface for momentum 
directions which are close to the normal of the interface. Indeed, expanding Eq. (16) in the vicinity of pj, 



un{\Al(p f )\, \Al(p f )\} = AI(P]) ~ 2 \A oc (p , j) (6p r V Pf ) A^p]) 



(29) 

and keeping only terms up to second order in 5pf, we get from (|16[), ([28|), ( |2S| ) the following range for the occurance 
of bound states: 



(dp r 



Aoo(py) 



Aoo(P^) 



< 



8 A 2 



A 2 (p n f ) 



(30) 



One can see from (g0() and condition \ip OD (p r j;) \ -C 1, that the right-hand side of Eq. (|3C|) is much smaller than 1 



B. Results for d-Wave Superconductors 

Now, we consider a tetragonal d x 2_ y a superconductor with a cylindrical Fermi surface whose z^-tods is parallel to 
the boundary plane. The order parameter is supposed to have the form 

A (pf, x) = A(x) cos(2<^ - 2a) , (31) 

where <fr is the azimuth angle between the Fermi momentum pf and the surface normal h (which is taken as the x 
direction). The angle a describes the orientation of the crystalline axis xq with respect to the normal of the boundary. 
For this case, Eq. fl2£|) for the bound state energy takes the form: 

-1- A 2 cos 2 (2a)- 



^Li cos 2 (2a) 



^ (Al cos 2 (2a) +4(1-2^ + 4^^) . (32) 

Consequently, the angular range for which bound states exist is limited by 

5^ 2 cos 2 (2a) 7T 
16tan(2a) ' °^ a ^4 ' (33) 

The quantity Aq = Ao((f> — 0) depends on the misorientation angle a as well as on the temperature T. According to 
(|33|), bound states for quasiparticles with momentum directions close to the normal of the boundary always appear, 
with the exception of two obvious crystalline orientations: a) a = where the order parameter is not suppressed at 
the interface {Aq vanishes for a — » 0), and b) a = ±45° where the order parameter for perpendicular incidence is 



According to (JT2[) , ( pi| ) the characteristic length describing the localization of bound states near the surface is 
of the order of v/ cos (^/(Acxj (Pf) sm (Voo (Pf))) ~ £CO cos 'Z 1 / sm (y'oo) cos(20 — 2a). For momentum directions, for 
which the bound levels lie close to the edge of the continuous spectrum, one gets from ( |l5|) | sin (tpoa(Pf)) 1^1- The 
characteristic length is in this case much greater than the temperature dependent coherence length £(T) (provided 
cosc/> is not very small). The factor cos(20 — 2a) in the denominator may result in an additional increase of the 
characteristic length. 

The results of numerical calculations for the angular dependence of the bound state energy £b(4>) are presented in 
Fig. 2 for various misorientation angles a. The temperature is T — 0.1 T c . The spatially varying gap amplitude A(x) 
was evaluated self-consistently (see Ref. [t]). Bound states with nonzero energy £b(4>) occur in a comparatively narrow 
angular region in the vicinity of the surface normal (Fig. 2a). These bound states are also described analytically 
above. Bound states which occur for lattice to surface orientations a > 37° (Fig. 2b) (it is sufficient to consider 
0° < a < 45°) correspond to values of the parameter \A 00 (pf)cos((p 00 (pf)) A (pf,(p oo (pf))\/ A max of the order of 



unity for which condition ( |20| ) is violated. They arise for trajectories for which A 00 (pf)A 00 (p^.) < 0. In general, there 

are also bound states for glancing angles, i.e., <f> w ±90°, but they give only a very small contribution to the tunneling 
current, and we do not plot them here. The limit for bound states with nonzero quasiparticle energy is determined 
by momentum directions for which discrete bound state levels merge in the continuum. Note that for misorientation 
angles 37° < a < 45° bound states appear for e = (midgap states), for e s» A^cf) = 0°) = A max cos(2a) 
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(perpendicular incidence) and for e w A oc (</) = 45°) = A max sin(2a) (0 = 45°). As a result one finds three peaks in 
the conductance (see inset of Fig. 6a). 

In Fig. 3 we display the dependence of the bound state energy for perpendicular incidence (4> = 0) on the misori- 
cntation angle a for T = 0. We included the function Aoo(0 = Q)/A max = cos(2a) to demonstrate that the bound 
states with nonzero energy lie very close to the band edge. 

The angle-resolved density of states i/(<f>, x — 0; e) at clean surfaces (p = 0) is shown in Fig. 4 for the misorientation 
angles a = 0°, 20° and 45°. For a (100) surface (a — 0°) (Fig. 4a) the order parameter is homogeneous up to the wall, 
and the surface density of states coincides with its bulk value. In particular, there are no surface bound states. The 
divergences presented in Fig. 4a lie at the edges of the continuum (e = |Aoo(0)|). For a lattice to surface orientation 
a = 20° (Fig. 4b) bound states occur near the surface. For perpendicular incidence (0 = 0°) there is a bound state 
with nonzero energy near the band edge, while for = 30 and 45° midgap states appear (cf. Refs. ||,[7|). In Fig. 
4c we display the surface density of states for an ideal (110) surface (a = 45°). For a quasiparticle moving along 
a trajectory with angle <f> — 45° we find both a midgap state and a bound state near the edge of the continuum. 
Both arise for an ideal boundary under the condition A 00 (pf)A 00 (p ) < 0. For this special surface orientation we 

have Aoo(p/) = — A OQ (p ) for all pairs of incoming and reflected trajectories. From Eqs. (13) and ( |l4| ) we find 
ip(jjf,x) = 7r — ip(p ,x), and at the surface ip(pf,0) — (p(pj.,Q) = n/2. On the other hand, we obtain from ( p"5| ) for 
bound states near the continuum <Poo(Pf) = 7r — <Poo(Pf) ~ or ir. Hence, the phase <p has to vary comparatively 
rapidly with spatial position x. For this reason, our analytical approach, based on the approximation (fL8|)-(p0|), breaks 
down for this special case. 

The influence of surface roughness is shown in Fig. 5 for a lattice to surface orientation a = 20°. We observe the 
broadening of the bound states with increasing roughness. The midgap state is smeared out faster than the nonzero- 
energy bound state (Figs. 5a, b). The point is that interference effects at the surface, which are of importance for 
the midgap states, are destroyed by diffuse scattering of quasiparticles at the boundary. On the other hand, for the 
existence of the nonzero-energy bound state only the suppression of the order parameter near the wall (even if it is 
very small) plays a role. For very rough surfaces (Fig. 5c) the midgap state is completely suppressed whereas the 
nonzero-energy bound state merges in the continuum. A similar behavior of midgap states at rough surfaces was 
observed in Ref. |[ 



IV. TEMPERATURE DEPENDENCE OF THE PEAKS IN THE DIFFERENTIAL CONDUCTANCE 

We consider a tunnel junction with transparency D <C 1. First, we focus on specularly reflecting interfaces between 
a normal metal and a clean singlet anisotropically paired superconductor. Let the normal of the junction barrier 
plane be directed along the x-axis (n || x). The superconductor occupies the left half-space (x < 0, index I) and the 
normal metal the right one (x > 0, index r). Then, the dissipative current in first order in the barrier transparency 
D may be represented byo 

I oo 

UV) = ~eA [ f^^D (p l f ) J "fe^^-tanli^^^Oje) . (34) 



Here, vg(p l pX — 0; e) is the angle resolved density of states (|J) of the superconductor on the left side of the tunneling 
barrier. It is calculated for zero transparency (D = 0). The voltage V = — < & r is the difference between the electric 
potentials of the left and right electrodes, and A is the area of the interface. Properties of the normal metal enter 
only implicitly via the function D(j)f), which is the probability of a quasiparticle with momentum pf to tunnel across 
the barrier. Note that the charge e is negative in our notation. 

Most experiments measure directly the differential conductance G = dl/dV. From Eq. (|34| ) we get 



G(V,T)-— J — D[p f ) J de cogh2 (£ _^ , (35) 



which reduces at zero temperature to 



2Ql v l . 

(36) 



G{V,T = Q) = 2eM J ^~^D (p> f ) u s (p l f , 0; eV) 
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Thus, tunneling experiments at low temperatures measure essentially the surface density of states integrated over the 
Fermi surface. As we will show below the surface density of states of anisotropic superconductors deviates substantially 
from its bulk value if the surface is pair breaking. Hence, tunneling experiments are extremely sensitive to special 
features of the surface (orientation and roughness). 

We notice, that in case of a real order parameter the density of states is symmetric with respect to the Fermi energy 
Ef, i.e., u(pf, R; e) — v(pf, -R; — e). Consequently, the conductance in Eq. (|35| ) is symmetric: G(V) = G(—V). Hence, 
the asymmetry of the conductance observed experimentally in Ref. |3^ is beyond the scope of this article. 



A. Surfaces without Pair Breaking 



Let us first consider surfaces which are not pair breaking, 
with its bulk value, 



Then, the density of states at the interface coincides 



v s (p/,0;e) 



y/e*-\A(p f )\- 



e(N-|A(p/)|) 



(37) 



At low temperatures (T <C \eV\) we get the following expression for the differential conductance (Ao(j>/) = A(pf, T 
0)): 



G « e A 



d 2 S l v l flX 
(2tt)3 |v^| 



D(p l f) 



de \eV\Q (jeVj - |A ( P/ )| + 2Tesgn(eVQ) 
cosh 2 (e) y/4eVTi+((eV)*-\Ao(jp f )\*) 



(38) 



In the case of an isotropic s-wave superconductor we obtain for T <C HeV^I — Ao| the well-known result for the 
conductance, 



G 



1 



\eV\ 



e(| e y|-A ) 



(39) 




At zero temperature the conductance diverges in the limit \eV\ — > A , which is directly associated with the singularity 
in the BCS density of states at e = Aq. It is worth noting that in the opposite limiting case, | |eT^| — Ao| <C T <C Ao, 
we obtain instead of ( |39| ) the following low-temperature anomaly in the conductance: 

(40) 

Hence, for the voltage \eV\ = Aq there is a divergence of the conductance, G oc 1/VT, in the low-temperature limit, 
due to the singularity in the BCS density of states. 

It has been known for a long time that anisotropy of the superconducting order parameter washes out structures 
in the integrated density of states, and correspondingly removes (or at least reduces) the singular behavior of the 
conductance (see Eq. (f36|)). This statement is valid, in general, in the absence of pair breaking at the surface. At 
the same time, in the presence of nodes of the order parameter on the Fermi surface, a nonzero subgap conductance 
appears. In particular, specific gapless and, in general, non-Ohmic behavior takes place in the limit of small voltages 
and low temperatures. For three-dimensional d-wave superconductor the singular behavior of the conductance vanishes 
in the absence of surface pair breaking. For the low -voltage behavior one has, depending on the crystal orientation of 
such a superconductor relative to the barrier plane:Q G oc |eV|/(i2jvAo) or G oc (|eV|/i?jvAo) ln(Ao/|eV|) under the 
condition T < \eV\ < A , and G oc {T/R N A Q ) or G oc {T/R N A ) ln(A /T) if \eV\ «T« A . 

Now, we consider a quasi two-dimensional Fermi surface, e.g., a Fermi cylinder with its principal axis parallel to the 
barrier plane. Let s = (sm, sj_) be a two-dimensional parameter set which describes the Fermi surface, i.e., Pf = Pf(s). 
As all quantities depend only on sii the order parameter can be approximated in the vicinity of a local maximum by 
= A max (l — bs?,). Substituting this expansion into (|3q ) we find the following singular contribution to the 
conductance (coming from the region bsf, < 1): 



G. 



ds± 



dp l f 



ds« 



dp\ 



Wf\ 



111 
1 



hi 



\eV\ 

T 
3^ 



T« \\eV\- A |, A , 
\eV\ - A | « T « A , 



(41) 
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where Ao = A max (T = 0). For the low voltage behavior one finds in this case in the presence of nodes of the order 
parameter G oc \eV\/R N A for T < \eV\ < A and G oc T/(R N A ) for \eV\ «T« A . 

So, in case of a d-wave superconductor with a cylindrical Fermi surface the low-temperature anomaly in the 
conductance as well as the singular dependence on the voltage turns out to reduce to a logarithmic singular- 
ity, G oc ln(T) , ln(| \eV\ — Ao|), which is much weaker compared to the s-wave case where we found G oc 

T- 1 / 2 , ({eV) 2 ~ AlY 1 ' 2 . 



B. Surfaces with Pair Breaking 



For ideally smooth surfaces with pair breaking the singular behavior of the conductance of an NIS junction with a 
d-wave superconductor may become comparable or even stronger than for junctions with isotropic s-wave supercon- 
ductors. The main reason for this are contributions from quasiparticle states bound to the barrier plane. 

If in the vicinity of a given voltage the dominating contribution to the conductance comes from bound states, one 
may consider only a singular part (a pole-like term) of the retarded Green function at the interface: 



9s (Pf, x = 



Qg (Pf) 



e - E B (pf) + irj 



(42) 



This yields a <5-peak in the angle resolved density of states 



■(Pf,x = 0;e) =Q g (p f )5{s-e B {Pf)) 



(43) 



Hence, the function eb (Pf ) describes the energy dispersion of the quasiparticles bound to the barrier. Note that Eb(pj) 
always has an extremum for the momentum direction along the normal of the boundary plane (pf\\h), provided the 
bound state occurs for this direction (see (p8j), Fig. 2 and Ref. [To]). Substituting ( |43| ) into (|35|) we obtain the following 
anomalous contribution to the conductance: 



G S (V,T) 



'IT 



d 2 S l v l - 



f,x 



(2irf 



DM) 



Qg [P 



cosh 



i(p' f )-eV 
IT 



(44) 



Consequently, the conductance rapidly increases with decreasing temperature, if the voltage is equal to the extremum 
of Eb(p f) ■ In case of a real order parameter bound states always appear in pairs with energies ±£b(pj) (see Eqs. 
(H), |J). Thus, we get from Eq. @) the relation G S {V) =G S {-V). 

Let us again consider a two-dimensional system. We assume that the dispersion of the bound states can be 
approximated around its maximum by £b( s ) — s B lax — fes| (b > 0). Then for \e g lax — \eV\ \ <Twe obtain from Eq. 
( [fif ) the following low-temperature anomaly in the conductance: 



G s 



(2V2-l) C (f)e 



dpK 



dpi 



(45) 



It is essential for the derivation of this result, that under the condition \s r . 



- \eV\\ < T only a region 6s 2 < T 

contributes to the integral over sn (due to the exponential decrease of the cosh _2 -function outside), and the condition 
T 6s 2 max permits to carry out the integration to infinity. Thus, making use of only the first two terms in the 
expansion of esiPf) near its maximum may be justified at sufficiently low temperatures. It is assumed for simplicity, 
that for all values 6s 2 < T the bound state still exists. All other quantities are supposed to change only slightly within 
the region &s 2 < T. 



Provided T < \e B lax 



\eV\\,bs- 



\,max 



we obtain instead of (|45j) the following dependence of the conductance on the 



voltage near the low-temperature peak: 



2e 2 AQ(e B lax - \eV\) f ds ± 
Vb ^/e^ ax -\eV\ J (2tt)- 



dp\ 



dsn 



dpi 



dsi 



r^D(P l f)Qg(P l f) 



(46) 



In the case of a three-dimensional Fermi surface we expect the function es(p/) to vary around its maximum as 



en Op/) 



b\S\ — bis\ (pi, 62 > 0). In contrast to the two-dimensional case eb(pj) now has an extremal point 
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instead of an extremal line. Then Eq. (Ejj) is valid only under the condition 62^2 max ^ ^i s i,moi ( e -S-; f° r a strongly 
elongated ellipsoidal Fermi surface) for temperatures \eg ax — \eV\\, b2s\ max <C T <C bis\ max . For lower temperatures 
r ^5 moi there are no temperature anomalies in the conductance. 

Due to the dispersionless character of midgap states = 0) their influence on the differential conductance 

differs from the influence of nonzero-energy states. We find from (H) 

v l , >0 

which is independent on the exact form and the dimensionality of the Fermi surface. The general analytical expression 

for the function Q g (p/J for midgap states via the spatially dependent order parameter as well as the explicit form of 

this function for a particular surface configuration are given in Ref. |l(]. In realistic systems the width of this zero-bias 
anomaly of the conductance is governed, apart from temperature, by the width of the broadened midgap delta-peak 
in the density of states. 



C. Numerical Results 



We have carried out numerical calculations of the conductance under various conditions. The angular dependence 
of the transparency is taken to be D(<f>) — Do cos 2 (0) where Dg is the transparency for perpendicular incidence and is 
related to the normal state resistance via Eq. (Q). The conductance G(V) at zero temperature is displayed in Fig. 6 for 
various misorientation angles a. Figure 6a shows the results for a smooth surface (p = 0), whereas in Figs. 6b, c the 
roughness parameter p is chosen to be 0.02 and 0.10, respectively. For a = there is no surface pair breaking and no 
quasiparticle states bound to the interface. In the case of an ideal boundary (Fig. 6a) the corresponding conductance 
has a logarithmic divergence (see Eq. (fi~L|)) at \eV\ — Aq from both sides of this point. In contrast to this peak, one 
can see asymmetrical peaks for a ^= 0. For finite misorientation angles a bound states, localized near the surface, 
arise. The conductance has a square root divergence at \eV\ = e r ^ ax (see Eq. (flq)) as well as a divergence at V — 
due to midgap states. Surface roughness (Figs. 6b, c) cuts these divergences and leads to finite peaks which become 
smaller and wider with increasing roughness. It is evident from these figures that the conductance peak at nonzero 
voltage does not measure the maximum gap A max but, to a good approximation, the (bulk) gap for perpendicular 
incidence A oc (p , j). We want to stress that the position of the peaks changes only slightly with increasing roughness. 
Although a direct comparison of our data with the results of Ref. is not possible (they consider the surface density 
of states and not the conductance) the influence of roughness is similar in both models. 

The low-temperature behavior of the conductance peaks is presented in Figs. 7a, b for smooth interfaces (p = 0) 
with misorientation angles a = 20° and 30°. The temperature dependence of the zero-bias anomaly G(0) is described 
in Fig. 7a, while in Fig. 7b the values of the peak at nonzero voltage G(V max ) are shown. We see, for example, that for 
a lattice to surface orientation a = 20° the zero-bias conductance G(0) is proportional with a good accuracy to 1/T 

up to T — 0.2T C , whereas the temperature dependence of the quantity G(V max ) is described by the fl/VT^-term 

only for a very small temperature range (T < 0.01 T c ). 

The influence of finite temperatures on the conductance is plotted in Fig. 8 (Fig. 9) for a lattice to surface orientation 
a = 20° (a = 45°). We present results for different degrees of interface roughness: (a) smooth interface (p = 0), (b) 
slightly rough interface (p — 0.02) and (c) higher degree of roughness (p — 0.10). The peaks are well pronounced only 
for T < 0.2 T c . Even at T = 0.05 T c (this corresponds to the temperature of liquid Helium if T c = 93 K is assumed) 
the reduction of the height of the conductance peaks is evident. At least at low temperatures, the effect of thermal 
smearing is more pronounced for the nonzero-bias peaks (Fig. 9a). Thus, we would like to encourage experimentalists 
to do tunneling measurements even below 4.2 K. 

Figure 10 compares the results for the conductance at T = using a constant order parameter (cf. Ref. |34|) 
with calculations including a self-consistent order parameter. In both cases (a = 20° and a = 45°) we find dramatic 
differences at finite voltages because of the appearance of nonzero-energy bound states due to the depletion of the order 
parameter near the interface. The presence of the zero-bias conductance peak is unaffected by the self-consistency of 
the order parameter, although disregarding surface pair-breaking results in overestimating of the weight of the peak 
up to 50% depending on the crystal orientation. 
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D. Summary 



In summary, we have studied in detail the appearance of surface bound states and their relevance for the tunneling 
spectrum of NIS junctions. We developed, for the first time, an analytical approach to the surface problem which 
allows an understanding of the presence of zero-energy as well as nonzero-energy bound states. Numerical calculations 
have been done for d-wave superconductors assuming a cylindrical Fermi surface. 

It is well known that zero-energy bound states are a robust phenomenon^ and will always occur whenever the real 
order parameter changes sign along a quasiparticle trajectory. SoJjar, it was not clear under which condition nonzero- 
energy bound states arise. They were only observed numerically.Bfl The result of our analytical consideration is that at 
least some of the nonzero-energy bound states are also of quite general (although not topological) origin and appear, 
for any (even tiny) suppression of the order parameter, within a narrow region around the surface normal (see Fig. 
2a). The reason is that the order parameter amplitudes for the incoming and the outgoing trajectories far away from 
the surface are approximately the same in this region, i.e., A oc (pf) « Aoo(p^). As a consequence, nonzero-energy 
states around the surface normal occur in the case of d-wave superconductors for all lattice to surface orientations 
with the exception of (100) surfaces (no suppression of the order parameter) and of (110) surfaces (A(pJ) = 0). In 
addition, we found numerically nonzero-energy bound states for misorientation angles between 37° and 45° as shown 
in Fig. 2b. These states are, however, beyond the validity of our analytical approach. 

Furthermore, we have studied the influence of bound states on the tunneling conductance, in particular at low 
temperatures (Figs. 6a). While zero-energy bound states always result in a 1/T divergence of the conductance for 
T — ► 0, the influence of nonzero-energy bound states depend on the geometry of the Fermi surface. For quasi 
two-dimensional Fermi surfaces the conductance diverges as 1/%/T at voltages determined by the extrema of the 
bound state energies. Within a few percent, these peaks coincide for a < 37° with the bulk order parameter for 
perpendicular incidence, i.e., \eV max \ s» |Aoo(p™)|. Only for higher misorientation angles (a > 37°) the conductance 
peaks correspond approximately to the bulk gap A max = max p/ |Aoo(p/)|. By contrast, low-temperature peaks in the 
conductance at nonzero-bias voltage are absent for three-dimensional systems. Our finite temperature calculations 
(Figs. 8a, 9a) demonstrate that the peaks can only be observed for temperatures below 0.2 T c . 

Finally, we have investigated the influence of surface roughness on the conductance (Figs. 6b, c, 8b, c, 9b, c). The 
peaks are broadened, but their positions remain essentially untouched. The broadening of the peaks with increasing 
surface roughness is more crucial for the zero-bias peaks. 

We conclude that a systematic experimental study of low-temperature peaks in the conductance of NIS junctions 
for a series of lattice to surface orientations would give valuable information on the anisotropy of the superconducting 
order parameter. We stress once more that nonzero-bias peaks give information on the gap amplitude in the bulk, 
even though the order parameter is suppressed at the surface. 
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FIG. 1. Schematic geometry of our interface model: a smooth interface covered by a dirty layer on each side. Also shown 
are the four Fermi momenta which are involved in the scattering process at the smooth interface (x = 0). 

FIG. 2. Angular dependence of the bound state energies eb{4>) near = 0° (a) and near <f> — 45° (b) for various 
misorientation angles a. 

FIG. 3. Bound state energy for perpendicular incidence £b(4> = 0) as a function of the misorientation angle a. For 
comparison, we included the position of the band edge, A ao (4> — 0)/A max = cos(2a). 

FIG. 4. Angle-resolved density of states at a smooth surface (p = 0) for various trajectory angles (f>. The lattice to surface 
orientation is taken to be a = 0° (a), 20° (b) and 45° (c). 



FIG. 5. Angle-resolved density of states at rough surfaces (a = 20°) for the same trajectory angles <j> as in Fig. 4b. The 
roughness parameter p is 0.02 (a), 0.10 (b) and 0.50 (c). 

FIG. 6. Zero-temperature conductance as a function of the bias voltage for various misorientation angles a. The surface 
roughness is p = (a), 0.02 (b) and 0.10 (c). The inset of (a) demonstrates the presence of three conductance peaks for a = 40°. 

FIG. 7. Temperature dependence of the zero-bias anomaly G(V = 0,T) (a) and the nonzero-bias anomaly G(V = V max ,T) 
(b) for two lattice to surface orientation angles a = 20° and 30° . 



FIG. 8. Conductance G(V) for various temperatures T. The lattice to surface orientation is taken to be a = 20°. The 
figures display the effect of increasing roughness: p = (a), 0.02 (b) and 0.10 (c). 



FIG. 9. The same as in Fig. 8 but for a = 45°. 



FIG. 10. Influence of the self-consistency of the order parameter on the zero-temperature conductance of a smooth interface 
(p = 0). We show the data for two misorientation angles a = 20° and 45°. 
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